Asymptotic L1-decay of solutions of the porous medium equation to self-similarity

We consider the flow of gas in an N -dimensional porous medium with initial density v0 (x) ≥ 0. The density v(x, t) then satisfies the nonlinear degenerate parabolic equation vt = ∆v m where m > 1 is a physical constant. Assuming that (1 + $2 )v0 (x) dx < ∞, we prove that v(x, t) behaves asymptotically, as t → ∞, like the Barenblatt-Pattle solution V ($, t). We prove that the L1 -distance decays at a rate t 1/((N+2)m−N) . Moreover, if N = 1, we obtain an explicit time decay for the L∞ distance at a suboptimal rate. The method we use is based on recent results we obtained for the Fokker-Planck equation [2], [3]

Convergence of the mass-transport steepest descent scheme for the sub-critical Patlak-Keller-Segel model

Variational steepest descent approximation schemes for the modified Patlak-Keller-Segel equation with a logarithmic interaction kernel in any dimension are considered. We prove the convergence of the suitably interpolated in time implicit Euler scheme, defined in terms of the Euclidean Wasserstein distance, associated to this equation for sub-critical masses. As a consequence, we recover the recent result about the global in time existence of weak-solutions to the modified Patlak-Keller-Segel equation for the logarithmic interaction kernel in any dimension in the sub-critical case. Moreover, we show how this method performs numerically in one dimension. In this particular case, this numerical scheme corresponds to a standard implicit Euler method for the pseudo-inverse of the cumulative distribution function. We demonstrate its capabilities to reproduce easily without the need of mesh-refinement the blow-up of solutions for super-critical masses

Particle interactions mediated by dynamical networks: assessment of macroscopic descriptions

We provide a numerical study of the macroscopic model of [3] derived from an agent-based model for a system of particles interacting through a dynamical network of links. Assuming that the network remodelling process is very fast, the macroscopic model takes the form of a single aggregation diffusion equation for the density of particles. The theoretical study of the macroscopic model gives precise criteria for the phase transitions of the steady states, and in the 1-dimensional case, we show numerically that the stationary solutions of the microscopic model undergo the same phase transitions and bifurcation types as the macroscopic model. In the 2-dimensional case, we show that the numerical simulations of the macroscopic model are in excellent agreement with the predicted theoretical values. This study provides a partial validation of the formal derivation of the macroscopic model from a microscopic formulation and shows that the former is a consistent approximation of an underlying particle dynamics, making it a powerful tool for the modelling of dynamical networks at a large scale

Convergence of an entropic semi-discretization for nonlinear Fokker-Planck equations in Rd

A nonlinear degenerate Fokker-Planck equation in the whole space is analyzed. The existence of solutions to the corresponding implicit Euler scheme is proved, and it is shown that the semi-discrete solution converges to a solution of the continuous problem. Furthermore, the discrete entropy decays monotonically in time and the solution to the continuous problem is unique. The nonlinearity is assumed to be of porous-medium type. For the (given) potential, either a less than quadratic growth condition at infinity is supposed or the initial datum is assumed to be compactly supported. The existence proof is based on regularization and maximum principle arguments. Upper bounds for the tail behavior in space at infinity are also derived in the at-most-quadratic growth case

Contractive probability metrics and asymptotic behavior of dissipative kinetic equations

The present notes are intended to present a detailed review of the existing results in dissipative kinetic theory which make use of the contraction properties of two main families of probability metrics: optimal mass transport and Fourier-based metrics. The first part of the notes is devoted to a self-consistent summary and presentation of the properties of both probability metrics, including new aspects on the relationships between them and other metrics of wide use in probability theory. These results are of independent interest with potential use in other contexts in Partial Differential Equations and Probability Theory. The second part of the notes makes a different presentation of the asymptotic behavior of Inelastic Maxwell Models than the one presented in the literature and it shows a new example of application: particle's bath heating. We show how starting from the contraction properties in probability metrics, one can deduce the existence, uniqueness and asymptotic stability in classical spaces. A global strategy with this aim is set up and applied in two dissipative models

Stability and asymptotic analysis of a fluid-particle interaction model

We are interested in coupled microscopic/macroscopic models describing the evolution of particles dispersed in a fluid. The system consists in a Vlasov-Fokker-Planck equation to describe the microscopic motion of the particles coupled to the Euler equations for a compressible fluid. We investigate dissipative quantities, equilibria and their stability properties and the role of external forces. We also study some asymptotic problems, their equilibria and stability and the derivation of macroscopic two-phase models

Refined asymptotics for the subcritical Keller-Segel system and related functional inequalities

We analyze the rate of convergence towards self-similarity for the subcritical Keller-Segel system in the radially symmetric two-dimensional case and in the corresponding one-dimensional case for logarithmic interaction. We measure convergence in Wasserstein distance. The rate of convergence towards self-similarity does not degenerate as we approach the critical case. As a byproduct, we obtain a proof of the logarithmic Hardy-Littlewood-Sobolev inequality in the one dimensional and radially symmetric two dimensional case based on optimal transport arguments. In addition we prove that the onedimensional equation is a contraction with respect to Fourier distance in the subcritical case

Infinite time aggregation for the critical Patlak-Keller-Segel model in R2

We analyze the two-dimensional parabolic-elliptic Patlak-Keller-Segel model in the whole Euclidean space R2. Under the hypotheses of integrable initial data with finite second moment and entropy, we first show local in time existence for any mass of "free-energy solutions", namely weak solutions with some free energy estimates. We also prove that the solution exists as long as the entropy is controlled from above. The main result of the paper is to show the global existence of free-energy solutions with initial data as before for the critical mass 8 Π/Χ. Actually, we prove that solutions blow-up as a delta dirac at the center of mass when t→∞ keeping constant their second moment at any time. Furthermore, all moments larger than 2 blow-up as t→∞ if initially bounded

A decision-making Fokker-Planck model in computational neuroscience

Minimal models for the explanation of decision-making in computational neuroscience are based on the analysis of the evolution for the average firing rates of two interacting neuron populations. While these models typically lead to multi-stable scenario for the basic derived dynamical systems, noise is an important feature of the model taking into account finite-size effects and robustness of the decisions. These stochastic dynamical systems can be analyzed by studying carefully their associated Fokker-Planck partial differential equation. In particular, we discuss the existence, positivity and uniqueness for the solution of the stationary equation, as well as for the time evolving problem. Moreover, we prove convergence of the solution to the the stationary state representing the probability distribution of finding the neuron families in each of the decision states characterized by their average firing rates. Finally, we propose a numerical scheme allowing for simulations performed on the Fokker-Planck equation which are in agreement with those obtained recently by a moment method applied to the stochastic differential system. Our approach leads to a more detailed analytical and numerical study of this decision-making model in computational neuroscience

Critical mass for a Patlak-Keller-Segel model with degenerate diffusion in higher dimensions

"Vegeu el resum a l'inici del document del fitxer adjunt"